Escape Probability and Mean Residence Time 
in Random Flows with Unsteady Drift * 

James R. Brannan, Jinqiao Duan and Vincent J. Ervin 

"a^': 

0> I Clemson University 

^_, ; Department of Mathematical Sciences 

Q ! Clemson, South Carolina 29634, USA. 
^r^ ■ E-mail: duan@math.clemson.edu 

^ : September 1, 1999 

o: 

O; 

^~~' ! Abstract 

^^ ' We investigate fluid transport in random velocity fields with un- 

p ■ steady drift. First, we propose to quantify fluid transport between 

flow regimes of different characteristic motion, by escape probability 
and mean residence time. We then develop numerical algorithms to 
solve for escape probability and mean residence time, which are de- 
,^ ', scribed by backward Fokker-Planck type partial differential equations. 

J-; I A few computational issues are also discussed. Finally, we apply these 

^ ■ ideas and numerical algorithms to a tidal flow model. 
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1 Introduction 

The Lagrangian view of fluid motion is particularly important in geophysical 
flows since only Lagrangian data can be obtained in many situations. It is 
essential to understand fluid particle trajectories in many fluid problems. 

Stochastic dynamical systems arise as models for fluid particle motion in 
geophysical flows with random velocity field 

X = f{x,y,t) + a{x,y)wi, (1) 

y = 9ix,y,t) + b{x,y)w2, (2) 

where wi{t),W2{t) are two real independent Brownian motion processes, f,g 
are deterministic drift part, and a, b are the intensity coefficients of diffusive 
noise part, of the velocity field. Note that the generalized derivative of a 
Brownian motion process is a mathematical model for "white noise". For 
general background in stochastic dynamical systems, see [Hiniini- 

Deterministic quantities, such as escape probability (from a fluid domain) 
and mean residence time (in a fluid domain), that characterize stochastic dy- 
namics can be computed by solving the Fokker-Planck type backward partial 
differential equations. 

In a previous paper 0, when the drift is steady, i.e, f,g do not depend 
on time, we have quantifled fluid transport between flow regimes of different 
characteristic motion by escape probability and mean residence time; devel- 
oped methods for computing escape probability and mean exit time; and 
applied these methods in the investigation of geophysical fluid dynamics. In 
this paper, we further consider the case of unsteady or nonautonomous drift 
f{x, y, t),g{x, y, t), develop a numerical algorithm for computing escape prob- 
ability and mean exit time, and demonstrate the application of this approach 
to a tidal flow model. 

2 Stochastic Dynamics 

of Fluid Particle Motion 

For a planar bounded domain D, we can consider the exit problem of ran- 
dom solution trajectories of (H))-© from D. To this end, let dD denote the 
boundary of D. The residence time of a particle initially at {x, y) at time t 
inside D is the time until the particle flrst hits dD (or escapes from D). The 



mean residence time r(x, y^t) = t + u{x, y, t), where u{x, y, t) satisfies [H] 

ut + -a^U:,:^ + -b'^Uyy + f{x,y,t)u^ + g{x,y,t)uy = -1, (3) 

u{x,y,t) = 0, ix,y,t)edDx{0,T), (4) 

uix,y,T) = 0, {x,y) e D, (5) 

where T is big enough, i.e., T > sup T{x,y,t) and here the sup is taken over 
all {x,y) in a compact set, i.e., the closure of the domain D. 

Note that u{x, y, t) = t{x, y,t) — t is the mean residence time after time 
instant t, or it quantifies how much longer a particle will stay inside D after 
we observe it at position {x, y) at the time instant t. 

Let r be a part of the boundary dD. The escape probability p{x,y,t) 
is the probability that the trajectory of a particle starting at position {x, y) 
and at instant t in D first hits dD (or escapes from D) at some point in V 
(prior to escape through dD — F), and p(x, y, t) satisfies [121 CH] 

1 1 

Pt + -a^Pxx + -h'^Pyy + f{x,y,t)p^ + g{x,y,t)py = 0, (6) 

pix,y,t) = 0, {x,y,t)e{dD-T)x{0,T), (7) 

pix,y,t) = 1, ix,y,t)ETx{0,T), (8) 

p{x,y,T) = 1, {x,y)eD. (9) 

Note that p{x, y, t) depends on F, and it may be better denoted as pr{x, y, t). 
Suppose that initial conditions (or initial particles) are uniformly dis- 
tributed over D. The average escape probability P{t) that a trajectory will 
leave D along the subboundary F at time t, before leaving the rest of the 
boundary, is given by [E] 

^(^) = ^ f f Pi^,y,t)dxdy, (10) 

\D\ J JD 

where l-Dl is the area of domain D. 



3 Numerical Approaches 

For the backward type of partial differential equation (jH)), we reverse the time 

s = T -t 



Then the mean residence time u{x,y,s) (we still use the same notation) 
satisfies 

_ 1 2 1,2 

Us 0, Uxx ~r D ^yy 

+ fix,y,T- s)ux + g{x,y,T- s)uy + l, (11) 

uix,y,s) = 0, ix,y,s)edDx{0,T), (12) 

u{x,y,0) = 0, {x,y)eD. (13) 

Similarly for the escape probability p{x, y, s) (we still use the same nota- 
tion), we have 

1 1 

Ps = ■^a^Pxx + -b^Pyy + fix,y,T- s)px + g{x,y,T- s)py,{U) 

p{x,y,s) = 0, {x,y,s)e{dD-T)x{0,T), (15) 

pix,y,s) = 1, {x,y,s)ETx{0,T), (16) 

p{x,y,0) = 1, {x,y)eD. (17) 

A piecewise linear, finite element approximation scheme ^ was used for 
the numerical solutions of the escape probability p{x,y,s), and the mean 
residence time u{x, y, s), described by the parabolic equations (jUj), and (fT^ . 
respectively. By transforming back to original time t = T—s, we get p{x, y, t), 
and u{x,y,t). We have used a few different time-discretization schemes, 
including the implicit backward in time, and Crank- Nicholson scheme j3]. 
The code works also for boundary defined by a collection of points lying 
on the boundary. A piecewise cubic splines were constructed to define such 
boundaries. 

4 Application to a Tidal Flow 

To demonstrate the above ideas and numerical algorithms, we consider a 
tidal flow model. This flow model is very idealistic and here we just use it 
as an illuminating example. 

Oscillatory tidal water motions dominate a large part of the coastal re- 
gions. Beerens and Zimmerman [2j considered a tidal flow model with veloc- 
ity field 

Uo = TTsm^Tix) cosi^Tvy) + 7v\cos{2iTt), (18) 

Vo = — 7rcos(7rx) sin(7r?y), (19) 



where A is a parameter measuring the intensity of the tidal wave 7rcos(27rt). 
We take < A < 3 as used by Beerens and Zimmerman [2j. 

As pointed out by Beerens and Zimmerman |2], it is essential to include 
more complicated temporal modes in this model in order to describe more 
realistic tidal flows. We include random temporal modes or white noise in 
this model, i.e., we consider a tidal flow model with unsteady drift part and 
random diffusive part: 

Uo = 7rsin(7rx) cos{Try) + 7rAcos(27rt) + V2ewi, (20) 

Vo = —ncos{7rx)sm{ny) + V2ew2, (21) 

where e > is the constant intensity of the white noise. We assume that the 
random temporal modes are weaker than the time-periodic mode and so we 
take < e < 0.1 in the following simulations. 

We study the transport of fluid particles in this tidal flow model. The 
equations of motion of fluid particles in this tidal flow are 

X = TT sin(7rx) cos(7r|/) + 7rAcos(27rt) + v2et(;i, (22) 

y = —TT cos{iTx) sm{iTy) + V2ew2, (23) 

where Wi{t),W2{t) are two independent Brownian motion processes. The 
unperturbed flow, with no temporal periodic "tidal" mode cos(27rt) and no 
temporal white noise modes wi,W2 of this tidal flow model, is the so-called 
cellular flow: 

X = 7rsin(7ra;) cos(7ry), (24) 

y = —TT cos^TTx) sm{ny) . (25) 

Figure H shows the phase portrait of this unperturbed flow P^ - (j^3jl . 

The partial differential equations for mean residence time u of fluid par- 
ticles in a fluid domain D, and for the escape probability p of fluid particles 
cross a subboundary F of D, are the following (in reversed time s = T — t), 
respectively, 

+ [n sin(7rx) cos{7Ty) + vrA cos(27rT — s)]ux 

— ncos^nx) sin^Try) Uy + 1, (26) 
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Figure 1: Unperturbed flow 



Ps = e(Pxx + Pyy) 

+ [n sin(7ra;) cos(7ry) + nX cos(27rT — s)]px 
— IT coslnx) sm{7ry)py . 



(27) 



In practical numerical simulations, the "final" time T should be taken big 
enough so that the solutions u, p do not change within a reasonable tolerance. 
To do so, we monitor the mean-square difference of u (and also p) at time 
T and T + 1. When this difference is within a reasonable tolerance (we use 
0.001), we take the T as the "final" time; otherwise, we increase the value of 
T and do the simulation again, until the tolerance criterion is met. 

We take a fiuid domain D to be a typical cell, i.e., the unit square, in the 
unperturbed fiow; see Figure El 

Unlike the stochastic systems with steady drift as studied by Brannan, 
Duan and Ervin ^.3j, the mean residence time and escape probability depend 
on time (although the change is small in this specific example); see Figures 
El in El and ini All these plots are for A = 1 and e = 0.1. The mean residence 
time and escape probability for < A < 3 and < e < 0.1 display similar 
features. 

In this fiow model, it turns out that the average escape probability cross- 
ing the top boundary y = 1 does not change much with time, with value 







Figure 2: A fluid domain D 




Figure 3: Mean residence time after time t = 10 








Figure 4: Mean residence time after time t = 50.5 








Figure 5: Escape probability crossing the top boundary y 
t = 10 



1 after time 








Figure 6: Escape probability crossing the top boundary y 

t = m 



1 after time 



around 0.2639. We observe similar features for crossing other three side 
boundaries. 

In the tidal flow model (fTH j) -(fT ^ with only time-periodic "tidal" mode 
7rAcos(27rt), Beerens and Zimmerman j2] found that there are "islands" in 
the tidal flow and the fluid particles trapped in such islands will never escape, 
for < A < 3 as used by Beerens and Zimmerman |[2j. This phenomenon 
is common in non-dissipative, Hamiltonian planar systems. Although the 
dissipation in the oceans is small, "no matter how small the dissipation is, 
the (oceanic) fluid has substantial time to experience the action of dissipative 
forces" J2]- So this "islands" phenomenon does not appear to be likely in 
realistic tidal flows; see more physical discussions in [3 El E] • 

While in our tidal flow model (j^n|) - ()2H) with both time-periodic "tidal" 
mode 7rAcos(27rt) and temporal white noise modes, all fluid particles will 
eventually escape from any fluid domain in finite time after we first observe 
them; see Figures El IH This feature is true for any e > 0. It appears 
that our stochastic tidal model is a little more realistic than Beerens and 
Zimmerman's model Pj. We remark again that we use this simple tidal flow 
model to demonstrate the applications of mean residence time and escape 



probability. 

In summary, in this paper we have discussed the quantification of fluid 
transport between flow regimes of different characteristic motion by escape 
probability and mean residence time, developed numerical algorithms to solve 
for escape probability and mean residence time, and applied these ideas and 
numerical algorithms to a tidal flow model. 
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